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ABSTRACT 

We present the results of one year of pulse timing observations of PSR J1740-5340, an eclipsing millisecond 
pulsar located in the globular cluster NGC 6397. We have obtained detailed orbital parameters and a precise 
position for the pulsar The radio pulsar signal shows frequent interactions with the atmosphere of the companion, 
and suffers significant and strongly variable delays and intensity variations over a wide range of orbital phases. 
These characteristics and the binary parameters indicate that the companion may be a bloated main-sequence star 
or the remnant (still filling its Roche lobe) of the star that spun up the pulsar In both cases, this would be the first 
binary millisecond pulsar system with such a companion. 

Subject headings: globular clusters: individual (NGC 6397) — pulsars: individual (PSR J1740-5340) — binaries: 
close 



1. INTRODUCTION 

The milli second pulsar J1740- 5340 in the globular cluster 
NGC 6397 ( D'Amico et al. 2001 ) is a member of a binary sys- 
tem with a relatively wide orbit of period 1.35 days, a compan- 
ion with mass > 0. 19 M0, and is eclipsed for about 40% of 
its orbit at a frequency of i' = 1.4 GHz. In much tighter eclips- 
ing bin aries, like those containing PSRs B 1957-1-20 (Fruchtere t 



al. 1990), B 1744- 24A (^.yne et al. 19901 |Ni ce & Thorsett 1992| ) 
and J205 1-0827 ( Stappers et al. 1996 ), all having companions 
of mass Mc < O.IOM0, the eclipses are believed to be caused 
by a wind resulting from ablation of the companion by a rel- 
ativistic pulsar wind. PSR J1740-5340 is in a wide orbit and 
the wind energy density impinging on the companion is signif- 
icantly less than that estimated for the close eclipsing systems, 
and is unlikely to drive a wind of sufficient density off a degen- 
erate companion. In the discovery paper, D'Amico et al. (2001) 
argued that the companion could be an unevolved star of mass 
comparable to the turnoff mass of the cluster, ^ 0.8 Mq, releas- 
ing a wind sufficiently dense to produce the observed eclipses. 
Such a hypothesis, which would make PSR J1740-5340 an un- 
usual system, can be ultimately tested with the optical identifi- 
cation of the pulsar companion. 

In this Letter we report on the first year of pulse timing obser- 
vations of this pulsar. We provide timing parameters, including 
the position, and we report on further observational evidence of 
the interaction of the pulsar signal with the companion atmo- 
sphere. 

2. OBSERVATIONS AND RESULTS 

Regular timing observations have been made since 2000 July 
with the Parkes radio telescope, using the center beam of the 
multibeam receiver at 1.4 GHz. The o bserving system is th e 
one used in the discovery observations (D'Amico et al. 2001). 



bank having 512 x 0.5 MHz frequency channels for each of 
two linear polarizations. After detection, the signals from in- 
dividual channels are added in polarization pairs, integrated, 1 
bit-digitized every 125 /is, and then recorded on magnetic tape. 
Pulse times of arrival (TOAs) are determined offline by fitting a 
standard pulse profile having a high signal-to-noise ratio to the 
observed profiles and analyzed usin g the program T EMPO* and 
the DE200 solar system ephemeris ( Standish 1982| ). A simple 
Keplerian model of a binary system was used in the analysis. 

For about six months we routinely collected data on this pul- 
sar several times a week using an integration time of 30 min for 
each observation. From the orbital parameters obtained dur- 
ing the confirming observations in 2000 July, we knew that the 
pulsar was eclipsed in the orbital phase interval 0.05-0.45, so 
we have generally made observations at epochs when the pre- 
dicted orbital phase was outside this interval. In this way we 
have obtained a reasonable number of good quality integrated 
pulse profiles (with signal-to-noise ratios in the range 10-25 
and typical timing errors in the range 10^0 /is), while some 
low quality profiles were rejected. Most of the latter were ob- 
tained at epochs relatively close to the predicted eclipse edge, 
but some were taken far from the nominal eclipse region, at 
orbital phases where the pulsar was expected to be detected 
clearly. We also noted that, while we could obtain a coher- 
ent timing solution over several months, the resulting post-fit 
rms scatter of the residuals was large compared to the typical 
nominal timing uncertainties of individual TOAs, suggesting 
the presence of systematic effects. 

In order to probe the hypothesis that this behavior could be 
due to the interaction of the pulsar radiation with the compan- 
ion atmosphere, we have recently initiated a series of long ob- 
servations starting close to the eclipse edge near the descending 
node, and lasting for a substantial fraction of the orbit in the 
non-eclipsing region. 



The effects of interstellar dispersion are minimized with a filter 

'http://pulsar.princeton.edu/tempo. 
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Fig. 1. — Observed signal intensity at 1.4 GHz as a function of orbital phase and pulsar phase for four long observations away from the nominal eclipse region. 
The data are processed in contiguous integrations of 120 s duration, (a) ~ 10 hr observation starting on 2000 December 3 at 22:23 UT; (b) ~ 4 hr observation starting 
on 2001 March 1 at 22:32 UT; (c) ~ 7 hr observation starting on 2001 March 5 at 18:3 1 UT; (d) ~ 1 1 hr observation starting on 2001 March 12 at 15:02 UT. 



The orbital period is ^ 32 hr and the pulsar is visible at 
Parkes for ^ 11 hr on a given day. Thus far we have obtained 
four relatively long observations. The results are shown in Fig- 
ure 1. In one case (Fig. la) the pulsar signal was relatively 
stable in intensity and pulse phase, while in the other three 
cases (Figs, lb, c and d) significant excess propagation delays 
of up to ^ 3 ms and strong intensity variations are observed. 
In particular, in the case of Figure Ic, the pulsar signal was 
always rather weak and the pulse was about two times wider 
than usual, covering about 30-50% of the rotational phase. The 
correlation between delay peaks and flux reduction suggests 
that interstellar scintillation is not responsible for this behav- 
ior. These fluctuations are comparable to those observed at 
a lower frequency (- 400 MHz) in PSR B1957H-20 near the 
eclipse edges around orbital phases 0.2 and 0.3 (F ruchter et 
al. 1990). In PSR J205 1-0827 ( [Stappers et al. 199(^ ) flux den- 
sity variations and variable propagation delays are seen at 1 .4 
GHz when the pulsar signal is detectable through the eclipse 
region. PSR B1744-24A displays similar irregularities also far 
from the usual eclipse r egion, with typical delays of 300/is 
(Nice & Thorsett 1992). Instead, much larger excess delays, 
^ 1 ms, are observed at nearly all orbital phases in the case of 
PSR J1740-5340. These observations suggest that the pulsar is 
in fact orbiting within an extended envelope of matter released 
from the companion. 

We investigated the frequency-dependent behavior of the ob- 
served delay events by splitting the 256-MHz bandwidth avail- 
able at 1.4 GHz into two adjacent subbands. As can be seen 
from Figure 1, significant pulse broadening and signal-to-noise 
reduction is usually associated with large delays. We have 
therefore limited our analysis to small delays only. Figure 2 
shows an example of the measured excess delays in the two 



subbands. The small panel also shows these measurements 
plotted against each other The best-fitting straight line cor- 
responds to At oc ;y-2 02±o.3{) j.jjg excess delays, consistent 
with the i''^ dependency expected if the delays are due to dis- 
persion in an ionized medium. If entirely due to dispersion, 
the excess propagation delays of up to ~ 3 ms visible in Fig- 
ure 1 would correspond to electron column density variations 
of - 1 .5 X 10"* Af_3 cm-2, or ADM - 0.47 Af_3 cm'^pc ( Af_3 
being the delay at 1.4 GHz in ms). This corresponds to a pulse 
broadening of 0.10 P Af_3 = 0.36 A f_3 ms over the receiver 
bandwidth (to be compared with an intrinsic pulse FWHM of 
0.17 P) and may be responsible for most of the signal attenua- 
tion and pulse broadening observed near the delay peaks. 

From Figure 1 we can also see that far from the eclipse re- 
gion the typical time scale of the rise to maximum delay is of 
order 1-2 hr. This indicates that observations of ^ 3 hr are re- 
quired in order to identify reliably a non-delayed portion of an 
observation and to obtain a bona fide TOA. This explains why 
a timing solution mainly based on relatively short integrations 
may be affected by systematic effects. 

In order to obtain a reliable phase-connected timing solution 
for this pulsar we have selected by visual inspection among all 
the data collected at 1.4 GHz in the last 12 months only those 
TOAs with no indication of delay and those corresponding to 
profiles not affected by pulse broadening. In this way we have 
obtained a solution with a post-fit rms residual of ~ 120 /is. 
The corresponding model parameters are reported in Table 1 . 
Figure 3a shows the residuals for all TOAs collected whereas 
Figure 3b shows those for the selected TOAs only, all relative 
to the model in Table 1. The dispersion measure (DM) was 
fitted using TOAs from two adjacent subbands from the sta- 
ble long integration shown in Figure la. A fit for orbital ec- 
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centricity may be biased because only ^ 60% of the orbit is 
sampled. Nevertheless a formal fit results in an upper limit of 
e < lO"'* (30-), and in deriving the parameters listed in Table 1 
we have fixed e = 0. Further observations may better constrain 
this limit. Similarly, a longer data span will be needed for in- 
vestigating possible changes in the orbital parameters, as seen 
e.g. in PSR B 1957+20 (Nice, Arzoumanian & Thorsett 2000) . 



3. DISCUSSION 

PSR J1740-5340 has the longest orbital period (32 hr) and 
the most massive minimum companion mass (0.19 M0, as- 
suming a neutron star mass of 1.40 Mq) among the 10 known 
eclipsing pulsars (Nice, Arzoumanian & Thorsett 2000; Camilo 
et al. 2000). When combined with the typical duration of the 
eclipse at 1.4 GHz (^13 hr), these parameters imply an eclipse 
radius Re larger than the orbital separation a for any value of 
the orbital inclination /. Re is defined as the chord subtended 
on the orbital circumference of radius a by the angle between 
the orbital phase of eclipse ingress (or egress) and the orbital 
phase 0.25. Table 2 shows the relevant parameters associated 
with different companion types: a 0.19 M© He-white dwarf 
(WD) (corresponding to / = 90°) and a range of main sequence 
(MS) stars spanning the mass interval M„s ^ 0.19-0.8 Mq 
(the latter value being the maximum compatible with the age of 
NGC 6397). As can be seen. Re 7.5 R© > a ^ 6.5 R©. 
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Fig. 2. — Excess group delays measured in two 128 MHz-wide bands cen- 
tered at 1454 MHz {filled symbols, connected by a solid line) and at 1326 MHz 
{open symbols, connected by a dotted line) for an event occurring at orbital 
phases 0.65—0.75 on 2001 March 1 (see Fig. lb). The delays in the two sub- 
bands are fitted (inserted panel) with a straight line of slope 1.21 ±0.04, corre- 
sponding to a frequency dependence for the delays (At oc u^"/ ) with a best-fit 
power-law index 7 = 2.02 ± 0.30. 

From the observation in Figure la, we infer that there are 
no systematic electron column density variations over a wide 
orbital phase interval. Hence the neutron star is not orbiting in- 
side a steady corona of ionized matter surrounding the compan- 
ion star. More likely, the pulsar is just skimming a large enve- 
lope of matter, whose clumpness and varying shape at distances 
from the companion comparable with Re could be triggered by 



the movement and the energy flux of the pulsar itself. The irreg- 
ular features appearing in the three observations of 2001 March 
(Figs. Ib-ld) could for example represent this kind of interac- 
tion. 

The value s of the Roche lobe radius Re span the interval 
1. 3-2.2 R© ( [Eggleton 1983| ). Because Re > a > Re, the mat- 
ter causing the eclipses escapes the gravitational influence of 
the companion star and therefore must be continuously replen- 
ished. The large position offset of PSR J1740-5340 with re- 
spect to the center of NGC 6397 (r ^ 0'.55) limits the un- 
known contribution to the observed period derivative P due 
to acceleration in the globular cluster potential. Assuming a 
King model (1962) for the central mass density, core radius 
rc = 0.'05 and ce ntral luminos ity density (in units of L©/pc^) 



log(pz,[0]) = 5.6 ( [Harris 1996[ ), distance d = 2.6 kpc (Reid & 
Gizis 1998), andamass-to-luminosityratioM/L = 3.3 (Pryor& 
Meylan 1993), we derive (see e.g. Camilo et al. 2000) |Pacc| < 
4TrG/{9c)(rJr)r,d(M/L)pL{0)P - lO'^". Even accounting for 
the uncertainties in some of the above parameters (e.g., rc), the 
observed P (Table 1) should be only slightly affected by Pace- 
Thus we can use P to infer the rotational energy loss of the pul- 
sar, £■ 1.4 X lO-'^ ergs s"'. A He-white dwarf of minimum 
mass M„,y/ = 0.19 M© with an effective temperatur e Th-j = 10'*K 



has radius R 



0.033 R© (iDriebe et al. 1998|, SerenelH et 



al. 2001). Because the radius decreases with increasing tem- 
perature, we use Riyd as an upper limit. At the distance of 
6.0 Rq, the power impinging on the degenerate companion 
would be Wwd ^ 10^" ergs s"'. This would allow release of 
an ionized wind from the degenerate companion at a maximum 
rate M„d,max = W^.d{R„d/GM„d) ^ 1 -6 x lO"'- Mgy"' . 
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Fig. 3. — (a) Timing residuals of all the TOAs collected (525 points derived 
from 46 observations of length ranging from ~ 30 min up to ~ 11 hr) relative 
to the timing model in Table 1. (b) Timing residuals of the 241 TOAs used to 
derive the phase-connected solution reported in Table 1 . The scales of the two 
panels are identical. The post-fit timing residuals in panel b, while smaller on 
average than those in panel a, ai'e still noticeably large compared to the nominal 
TOA uncertainties. This effect may be caused by small-scale electron density 
fluctuations not modeled in our present fits. 

Assuming isotropic emission, mass continuity implies that 
the companion star releases mass at a rate = 47ri?|p(/?£)v/. 
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Using the values in Table 2 and the density at the eclipse ra- 
dius inferred from the observation of the excess delays assum- 
ing completely ionized matter {p[Re] ^3 x 10~'**Af-3 g cm"-'), 
we obtain Mj. ^ 1.5 x 10~" Af_3 vy.g MQy"' where vy.g is the 
wind velocity at Re in units of 10^ cm s"' (typical order of 
magnitude of the escape velocity from the surface of the com- 
panion). Thus, only in the case of low terminal velocity of the 
wind (vf <^ 100 km s"', corresponding to improbable electron 
temperatures of <^ 200 K at Re) could the pulsar flux sustain the 
ablation. More massive WDs would further exacerbate these 
difficulties. 

The hypothesis of a MS companion appears more viable: 
such a companion could have been acquired as a result of an 
exchange interaction in the cluster core (Sigurdsson & Phinney 
1993; Heggie, Hut & McMillian 1996), a scenario possibly sup- 
ported also by the position of PSR J1740-5340 ^11 core radii 
from the center of NGC 6397. This object is the second farthest 
pulsar from a cluster core (in units of core radii) among the 36 
pulsars in globular clusters with precisely determined positions 
(after PSR B2127H-1 IC in M15, located at - 13 core radii). 

The unperturbed radius R^, of a MS star of mass M,„5 ~ 
0.19-0.8 Mq is ^ 6-30 times greater than the maximum pos- 
sible radius for a WD (Table 2). The energy required for trans- 
porting a unit mass from the surface of the companion to its 
Roche lobe radius Ri ^, is t/* ~ GM^{R~J -Rl]j- Given this, the 
corresponding maximum ablation rate for a MS companion is 

Mms,max ~ (Rms / Rwdf(U,,d /U,„,)M,,d.max ^ 200-6000M ,,d .max = 

(3-90) X 10-'° Mgy-i. Therefore, only 0.2-5% of the imping- 
ing power would be needed to sustain the inferred mass loss 
rate Mc- The bulk of the pulsar energy would thus be available 
for heating the companion photosphere and perhaps inflating its 
radius up to Roche lobe, in turn making the release of the wind 
even easier This model should be more easily applicable to a 
light MS star, as large convective envelopes favor bloating of 



the star (Podsiadlowski 1991 ). 



Another hypothesis for explaining the large mass loss rate 
Mc is that the companion is an evolved star that spun up the 
millisecond pulsar and that is presently undergoing the pre- 
sumed final stages of mass loss. In this case, mass accre- 
tion and spinning up of the pulsar would now be inhibited by 
the pulsar wind flux which could expel the matter overflowing 
from the Roche lobe of the companion (see for instance Bur- 
den et al. 2001 and references therein). If substantiated, this 
would be the first confirmed example of a recently born mil- 
lisecond pulsar. The apparent rather young spin-down age of 
PSR J1740-5340 (r,. = P/2P - 350 Myr, assuming as men- 
tioned earlier that Pace is negligible), and its surface magnetic 
field (B = 3.2 X ^'''[PP]'/^ - 8 x 10^ G) would support this 
model. However, the future evolution of the binary (companion 
evolving into a WD?) remains unclear. Detailed simulations of 
the system are required for comparing this or other more ex- 
otic scenarios (Burderi, D'Antona & Burgay 2001, in prepara- 
tion) with the parameters derived from radio and optical obser- 
vations. 

Ferraro et al. (2001) report the possible identification of the 
optical companion to PSR J1740-5340 with a bright variable 
star and Grindlay et al. (2001) have detected with Chandra the 
X-ray counterpart of the millisecond pulsar The availability of 
data in a large spectral band, from radio to optical and X-rays, 
should help understanding this exotic system and the mecha- 
nism(s) responsible for the radio eclipses. 

The Parkes radio telescope is part of the Australia Telescope which 
is funded by the Commonwealth of Australia for operation as a Na- 
tional Facility. We thank the staff at Parkes for their support of this 
project, and L. Burderi for carefully reading the manuscript and for 
stimulating discussions. N.D'A. and A.P are supported by the Min- 
istero della Ricerca Scientifica e Tecnologica (MURST). F.C. is sup- 
ported by NASA grant NAG5-9095. 
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Table 1 

Parameters of the PSR J1740-5340 system 



Parameter 




Value 


Right Ascension (J2()(«)) 




it 40"" 44?589(4) 
-53° 40'40"9(1) 


Declination (J 2000) 




Galactic longitude, / 




338f 164 


Galactic latitude, b 




-11° 96 


Dispersion Measure, DM (cm" 


"'pc) 


71.8(2) 


Period, P (s) 




0.00365032889720(1) 


Period derivative, P (10~^°) 




16.8(7) 


Period epoch (MJD) 




51917.0 


Orbital period (d) 




1.35405939(5) 


Projected semimajor axis (s) 




1.65284(7) 


Epoch of ascending node, Tasc 


(MJD) 


51749.710822(6) 


Eccentricity'', e 




< 10-^ 


Data span (MJD) 




51719-52116 


Post-fit rms (lis) 




114 


Mass function (Mq) 




0.0026442(3) 


Companion mass, (Mq) 




>0.19 



Note. — Numbers in parentheses represent three times the formal la uncertainties 

given by TEMPO in the least-significant digits reported. 

^AU other parameters have been derived setting e = (see text). 



Table 2 

Parameters for possible companion 



Parameter 


He-WD 


MS 


Companion mass, (Mq) 
Orbital inclination, i (° ) 


0.19 


0.19-0.8 


90 


90-17 


Orbital separation, a (Rq) 


6.0 


6.0-6.7 


Roche lobe radius, Rl (R0) 


1.3 


1.3-2.2 


Eclipse radius, R^ (Rq) 


7.1 


7.1-7.9 


Stellar radius, R(Rq) 


'-O.OS 


0.2-0.96 



